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An effective method to estimate composition amplitude of spinodal 
decomposition for atom probe tomography validated by phase field 
simulations 

Wei Xiong, John Agren, 1 and Jing Zhou 1 

Department of Materials Science and Engineering, KTH Royal Institute of Technology, SE-100 44 Stockholm, 
Sweden 

Reasonable evaluation of composition amplitude in spinodal decomposition is a challenge to microanalysis of 
atom probe tomography, especially at early stages when phase separation is not prominent. This impedes 
quantitative analysis of spinodal structure in atom probe tomography as well as comparison with simulated 
results from phase field simulations. We hereby report an effective method to estimate the composition ampli- 
tude by constructing an amplitude density spectrum. This method can sensitively determine the composition 
amplitude at early stages. In particular, it substantially bridges experimental and simulation techniques com- 
prising both discrete and continuum data in the study of spinodal decomposition. Moreover, it was found that 
the commonly adopted Langer-Bar-on-Miller method for atom probe analysis underestimates the composition 
amplitude of spinodal decomposition. Case studies have been performed on the Fe-Cr binary alloys. 
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I. INTRODUCTION 

In the study of phase separation, atom probe tomogra- 
phy (APT) is a preferential technique to determine segre- 
gation of solute atoms due to its atomic-scale resolution 
of the 3-dimensional (3D) distribution of atomsi~~— . In 
order to elucidate mechanism of phase transformation in 
materials, Monte Carlo (MC) technique is usually consid- 
ered as a simulation tool to provide data for comparison 
with APT, since both techniques generate the discrete in- 
formation of the solute atoms as the direct outputs. How- 
ever, MC is limited by its time and length scale, and thus 
difficult to predict microstructure evolution with the ex- 
perimental time scale. Alternatively, phase field method 
(PFM) is often employed, which is based on a contin- 
uum approximation developed by Hillert4 and Cahn and 
Hilliard 5 . Quantitative comparison for spinodal struc- 
ture between APT and PFM depends on reasonable es- 
timation of the composition amplitude and wavelength. 
Additionally, one may easily envisage that estimation 
of composition amplitude is more important than wave- 
length at the early stages of spinodal decomposition, 
when spinodal decomposition starts to form and wave- 
length is too short to determine precisely. Since mechan- 
ical failure of materials, e.g. embrittlement of stainless 
steels, usually occurs at early stages of spinodal decom- 
position, estimation of composition amplitude is of the 
highest interest in engineering applications 6,7 . Unfor- 
tunately, such a comparison is greatly impeded, since 
there is no method available yet to effectively analyze 
both continuum data from PFM and discrete data from 
APT. Therefore, a robust method to evaluate the com- 
position amplitude in both simulations and experiments 
is urgently needed. 

It is noteworthy that the most commonly adopted 
method for estimation of the composition amplitude 



in APT is the Langer-Bar-on-Miller (LBM) method^, 
which can be also applied to the discrete data from MC 
simulations^. In the LBM method, the estimation of the 
composition amplitude is based on a frequency distribu- 
tion histogram as shown in Figure 1(a) which stands for 
the probability of finding a certain composition in the 
analyzed volume of the alloy. In order to apply the LBM 
method, local composition should be estimated by sam- 
pling a certain number of atoms, a block, to calculate 
the average composition in that sampling block from the 
atom map as illustrated in Figure 1(a). In the LBM 
method, the frequency distribution histogram is repre- 
sented by a continuous distribution function p(x) given 
by a sum of two Gaussian distributions with the same 
composition variance a: 
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where x\ and x 2 are the compositions for the two Gaus- 
sian peaks, xo is the average composition in the analyzed 
volume of the sample. The experimental histogram thus 
is represented by 3 parameters x±, x 2 and <r, and the 
composition amplitude is (x 2 —x\). However, the tails of 
the fitted Gaussian curves can sometimes extend beyond 
the physical limits of and 1 at the composition axis 
(see Figure 1(a)). In fact, this may happen quite often 
when the sampling block for constructing concentration 
field is either comparable with or less than the half of 
the wavelength. Therefore, it artificially brings uncer- 
tainties into the analysis, and the fitted Gaussian distri- 
bution functions in the LBM method lose all statistical 
meaning. More importantly, when the tails extend to the 
unphysical region with the composition value below or 
above 1, there is no guarantee that Eq. (pQ) is consistent 
with the average composition of the alloy which would 
require that the integral from to 1 must be unity. It 
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FIG. 1. (a) The frequency distribution histogram of the sam- 
ple volume with 33.1 at.% Cr isothermally aged at 773 K for 
100 hours, constructed by using the block of 25 ions from 
the Cr atom map, which is partially shown on the right; (b) 
Constructed concentration field from the APT data. 



has been suggested that some manipulations, like rescal- 
ing and constraining the data over the whole composition 
range, may be considered in order to eliminate that the 
tail exceeds the composition limit o 10 ^ 1 . Obviously, the 
accuracy of the LBM analysis cannot be guaranteed by 
taking such an expedient. 

One shall notice that the LBM method absolutely loses 
its applicability to the continuum data. For instance, the 
frequency distribution histograms that can be generated 
from PFM simulations differ quite much from the ones 
measured by APT, see Figure 2(a) to (c). Since the phase 
field model is based on a continuum model and the com- 
position on each node of the simulation grid is obtained 
by solving the Cahn-Hilliard equation*^, the obtained 
frequency distribution histogram for a real system is un- 
likely to show any value at or 1 on the composition 
axis. Moreover, the frequency distribution histograms 
from phase field simulations usually exhibit much sharper 
changes than those obtained experimentally by APT with 
a bell curve shape shown in Figure 1(a)), and thus usu- 
ally cannot be fitted by the Gaussian distribution func- 
tions of the LBM method. A histogram more similar to 
the experimental one may be obtained by considering the 
simulated value in a node as an average in a block of a 
given size and assume that the atoms within that block 
are randomly distributed which leads to a binomial dis- 
tribution which may be approximated by a Gaussian if 

I 



the block is not too small. 

Accordingly, there is a strong need to find a versatile 
method for evaluating the composition amplitude of spin- 
odal decomposition for different techniques with discrete 
(APT or MC simulation) and continuum (PFM) infor- 
mation as the output. In this work, we propose such a 
new method, called Amplitude Density Spectrum (ADS), 
to estimate the composition amplitude in order to bridge 
different experimental and simulation tools for quantita- 
tive investigation of the spinodal microstructure evolu- 
tion. Furthermore, the case study on the Fe-Cr alloy will 
indicate advantages of the ADS method in this work over 
the LBM method. 

II. MATERIALS, EXPERIMENTS AND SIMULATIONS 

A binary Fe-Cr alloy was considered for testing the new 
method in the current work because of its importance to 
the investigation of spinodal decomposition in stainless 
steels 643 16 . Both APT and PFM were employed in the 
investigation. 



A. Atom probe tomography 

The Fe-Cr binary alloy with 37.8 at.% Cr was pro- 
duced by vacuum arc melting, and then homogenized 
at 1373 K for 2 hours under pure argon followed by 
quenching in brine. Afterwards, samples were annealed 
at 773 K for 100 hours, and then shaped into needle-like 
specimens for the standard two-stage electro-polishing. 
The analysis was performed under a local electrode atom 
probe (LEAP 3000X HRTM, Imago Scientific Instru- 
ments, USA) equipped with a reflectron for improved 
mass resolution. The ion detection efficiency is about 
37 %, and the experiments were made in voltage pulse 
mode (20 % pulse fraction, 200 kHz, evaporation rate 
1.5 %) with a specimen temperature of 55 K. It should 
be noted that the analysis of the concentration in the 
analyzed volume in APT was found to be 33.1 at.% Cr. 



B. Phase field simulation 

in order to verify the ADS method, a 3D phase field 
simulation on the Fe-35 at.% Cr alloy was performed 
based on the Cahn-Hilliard equation: 



dxcr 
dt 



V • X Cr XFe (x Cr M Fe + X Fe M Cr ) V 



(da 



\dxcr 



K • V XCr 
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The free energy G m is obtained from the most recent 
thermodynamic modeling by Xiong et a h 17 i 18 , while the 
parameters for describing atomic mobilities, Mc r and 
Mp ei are taken from Jonsson 19 . The simulation box 



is defined using uniform grids, which is equivalent to 
14 x 14 x 14 nm 3 . The composition gradient coefficient, 
ft, is defined as a function of the interatomic distance 
d (0.25nm), and regular solution parameter £7(12225.15 
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FIG. 2. Frequency distribution histogram and concentration 
field simulated by PFM at dimensionless time of (a) 2, (b) 6, 
and (c) 30. (d) Comparison of the analyzed results for com- 
position amplitude and its upper and lower limit at different 
simulation times using DFD and ADS methods. 

J/mol at 773 K according to the improved thermody- 
namic modeling 18 ), which equals to fld 2 /2. The semi- 
implicit Fourier-spectral method 12 is applied for simu- 
lation in a dimensionless form under periodic boundary 
conditions. 



III. RESULTS AND DISCUSSION 

A. Analysis of APT data using the LBM and ADS 
methods 

Firstly, the LBM method was applied to estimate the 
composition amplitude for the prepared Fe-Cr sample. 
Since the original information from APT is the position 
of individual atoms as shown in Figure 1(a), the local 
composition is estimated at first in a certain size of block, 
which was chosen to be 25 ions in this work. As a result, 
the frequency distribution histogram for a given bin size, 
see Figure 1(a), can be achieved by statistically count- 
ing the numbers of concentration values within a certain 
composition range, i.e. composition limit of each bin in 
the frequency distribution histogram. Alternatively, a 
concentration field was constructed by using a sampling 
block with a volume of 0.83 x 0.83 x 0.83 nm 3 , if a cu- 
bic shape is considered for the sampling block with 25 
ions. It is noteworthy that the frequency distribution 
histogram obtained in this way is almost identical to the 
previous one by counting the number of nodes of the 
concentration field with different concentration values in 
different concentration ranges. Using the LBM method, 
the frequency distribution histogram in Figure 1(a) can 
be fitted by two Gaussian distribution functions with the 
peaks at Cr content (mole fraction) of 0.249 and 0.416, 
respectively. As pointed out before, the frequency distri- 
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FIG. 3. (a) The ADS pattern constructed for the APT data 
of the sample with 33.1 at.% Cr. The analyzed composition 
limits by the ADS method are denoted as the red chained 
lines, while the ones from the LBM method are indicated as 
the green dashed lines; (b) the magnified part of the ADS 
pattern. The dashed lines indicate the average composition 
limits (0.243 and 0.418) before the first refinement of the 
ADS method, while the composition limits pointed by arrows 
(0.199 and 0.466) are the ones after the first refinement of 
the ADS method. The blue dots marked in the dashed circles 
should be eliminated before the first refinement. 



bution histogram in Figure 1(a) has a tail exceeding the 
composition limit to negative values, which indicates that 
the Guassian function cannot represent the frequency dis- 
tribution histogram satisfactory. Besides, one shall keep 
in mind that the outline of the histogram largely depends 
on its bin size. As a consequence, the choice of bin size 
is also a source of uncertainty in the LBM method. 

Based on the concentration field evaluated from APT, 
an ADS pattern could be directly constructed in Figure 
3. In fact, the ADS pattern is a ID plot of concentration 
values at each node in the concentration field following 
a spatial sequence over the whole analyzed volume. The 
overall ADS pattern normally contains much data points, 
and thus too dense to show clearly without large magni- 
fication, which can be found easily by comparing Figures 
3(a) and (b). Therefore, it is necessary to establish a 
criterion to select peaks and troughs for evaluating the 
composition amplitude. 

In this work, the evaluation of the composition am- 
plitude is based on the following empirical steps, and 
eventually leads to results which can be verified further 
by PFM. The method seems to be robust and capable of 
determining small amplitudes. First, all of the peaks and 
troughs are selected to obtain the average concentration 
for the peaks and troughs. As indicated in Figure 3(b), 
the average Cr concentration for peaks is 0.418, while the 
one for troughs is 0.243, according to the ADS data in 
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Figure 3(a). In order to achieve reasonable composition 
amplitude, some refinements are necessary. However, be- 
fore refinement, we shall remove those peaks and troughs, 
as the ones marked in Figure 3(b) with blue dashed cir- 
cles. The eliminated peaks are the ones lower than the 
average sample composition (0.331), and the excluded 
troughs are higher than the average. As the first step 
of the refinement, we only keep the peaks with the com- 
position value higher than the above computed average 
peak composition (0.418), and the troughs lower than 
the computed average trough composition (0.243). As a 
consequence, the first refinement generates two average 
concentration values as the updated criteria for the sec- 
ond refinement: 0.466 for peaks and 0.199 for troughs as 
indicated in Figure 3(b). Based on these two new crite- 
ria, further refinements will be performed iteratively as 
in the first refinement. In the present case, the average 
composition values from each refinement for upper and 
lower limits are denoted as the red solid squares in Fig- 
ure 4, respectively. Meanwhile, numbers of nodes at each 
step for refinement are plotted as the blue open circles. 

Furthermore, we found that the blue curves presenting 
the node numbers in terms of the refinement time are 
always varying smoothly. On the contrary, it is less com- 
mon to get a smooth profile for the curve representing 
average Cr concentration. Analogous to defining deflec- 
tions from a given baseline, e.g. onsets of transforma- 
tions in a signal curve from differential thermal analysis, 
a criterion is set to determine the composition for upper 
and lower limits of the composition amplitude, which is 
simply making a tangent cross (see black dot A in Fig- 
ure 4) of the blue curve for the number of nodes starting 
from the two ends as illustrated using the dashed lines 
in Figure 4. Afterwards, the composition limits can be 
read from the red star symbol B indicated on the com- 
position curve, which has the same value of the vertical 
coordinate as point A. It is suggested that, in order to 
minimize the artificial manipulation in data analysis, the 
curve of the number of nodes will not be fitted by a poly- 
nomial to obtain a tangent cross. Instead, the tangent is 
directly obtained from the linear extension from the first 
two ending points in the blue curve for the number of 
nodes as illustrated in Figure 4. 

As a result, the estimated composition amplitude of 
the sample volume with 33.1 at.% Cr annealed at 773 K 
for 100 hours is 0.348, with the upper composition limit 
of 0.507 and lower limit of 0.159 as indicated in Figure 
3(a). Apparently, the estimated composition amplitude 
from the ADS method (0.348) is much larger than the 
one by LBM (0.167). However, we shall notice that the 
average compositions of the peaks and troughs before 
refinements are 0.418 and 0.243 (see Figure 3(b)), which 
even show a larger composition difference, 0.175, than the 
evaluated composition amplitude from the LBM method. 
Obviously, large amount of gradient information has been 
included when computing the above two average compo- 
sition values. Hence, the composition amplitude should 
be at least larger than 0.175. This also means that the 
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FIG. 4. The upper (a) and lower (b) limits of the composition 
amplitude analyzed by the ADS method. Blue open circles 
indicate number of nodes at each refinement step, while red 
solid squares are for calculated composition limits in each re- 
finement. Points A and B have the same value on the vertical 
coordinate. 



LBM method has enormously underestimated the com- 
position amplitude. 

Moreover, one may realize that the size of the sam- 
pling block will affect both the frequency distribution 
histogram and the concentration field based on the dis- 
crete data from APT. A reasonable sampling size should 
well represent the periodicity of the connected structure. 
However, the manipulation is difficult for 3D cases, and 
deserves further discussion in another work. On the con- 
trary, the continuum data from phase field simulation 
can circumvent such a restriction due to sampling size. 
Therefore, the analysis of the phase field simulation re- 
sults is performed. 



IV. ANALYSIS OF PHASE FIELD SIMULATION 
RESULTS 

It should be emphasized that the analysis of the phase 
field simulation result is desired not only to approve the 
validity of the now suggested ADS method but also to 
show its advantages of analyzing the continuum data over 
the LBM method. 

In the present work, a typical feature of the frequency 
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distribution histogram from PFM is shown in Figure 2. 
In the early stages, as pointed out before, the outline of 
the frequency is similar to the one read from APT with 
two peaks overlapping shown in Figure 2(a). Figure 2(b) 
indicates a phase separation close to the medium stage. 
The peak on the Fe-rich side is pronounced, while the 
other peak is difficult to identify. However, when ap- 
proaching to the later stages as shown in Figure 2(c), the 
peak on the Cr-rich side starts to be distinguishable. So 
far, the only method to define composition limits of am- 
plitude is to determine the peak positions in frequency 
distribution histogram. If they are distinguishable as in 
the late stages (see Figure 2(c)), the composition ampli- 
tude will be the difference of the two composition lim- 
its corresponding to the peak position in the frequency 
distribution histogram. We named this as the direct fre- 
quency diagram (DFD) method. 

Accordingly, the analyzed results by both DFD and 
ADS methods are summarized in Figure 2(d). The 
dashed lines marked with the odd numbers correspond to 
the DFD analysis once the peaks are pronounced enough 
for getting the composition limit of amplitude. Over- 
all, the maximum difference of the composition ampli- 
tude between two different methods is 4.25 at.% Cr at 
dimensionless time of 7, when the peak of frequency dis- 
tribution histogram on the Cr-rich side starts to be dis- 
tinguishable. Apparently, the overall agreement between 
these two methods is fairly good. It should be noted that 
the definition of the composition amplitude in spinodal 
decomposition for real materials is blurred and depends 
on the determination method itself. An effective way 
to determine the composition amplitude should at least 
be capable to correctly and systemically demonstrate the 
variation of the composition amplitude along with the mi- 
crostructure evolution. According to the analysis by the 
ADS method, the tendency of the variation for composi- 
tion amplitudes is consistent with the one using the DFD 
analysis. Moreover, since the initial random noise set in 
the phase field simulation is 5.0 at.% Cr, it is encouraging 
to point out that the analysis based on the ADS methods 
obtains a value of 4.01 at.% Cr. Interestingly, the com- 
position amplitude in the phase field simulation shows a 
minimum close to the dimensionless time of 1, which is 
consistent with the spinodal theory. Because the initial 
input state only includes the random noise and needs a 
certain incubation time to form the connected structure 
with certain wavelength under the thermodynamic and 
kinetic conditions^. The above analysis with PFM has 
successfully demonstrated the validity of the ADS meth- 
ods. 

The continuum PFM has various advantages in sim- 
ulating the microstructural evolution. In order to inte- 
grate it with other experimental methods (e.g. APT) 
and simulation tools (e.g. MC, Molecular Dynamics), it 
is quite important to use a method, like ADS, to ana- 
lyze the results containing both discrete and continuum 
information in the related studies. 



V. CONCLUSIONS 

In summary, a new method called ADS is developed in 
order to evaluate the composition amplitude from both 
discrete and continuum information. According to the 
analysis on the 3D phase field simulation results, the 
ADS method is accurate and effective. Based on the 
analysis of the APT data, the ADS method is approved 
being more accurate than the LBM method when applied 
to the early stages of decomposition in Fe-Cr alloys, at 
which the mechanical properties may drastically change, 
and thus needs sensitive measurement to determine the 
intensity of the spinodal decomposition. As a tool of mi- 
croanalysis, the ADS method substantially bridges differ- 
ent experimental and simulation techniques for studying 
spinodal decomposition. 
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